Amendments to the Specification : 

Please replace paragraph [0005] as follows: 

[0005] Techniques for modelling fractured porous media are described in French 
Patents 2,575,947 2,757,947 and 2,757,957 filed by the assignee. The technique of the 
first patent relates to determination of the equivalent fracture permeability of a fracture 
network in an underground multilayer medium from a known representation of the 
network, allowing systematic connection of fractured reservoir characterization models 
to double-porosity simulators in order to obtain more realistic modelling of a fractured 
underground geologic structure. The technique of the second patent relates to 
simplified modelling of a porous heterogeneous geologic medium (such as a reservoir 
crossed through by an irregular network of fractures for example) in the form of a 
transposed or equivalent medium so that the transposed medium is equivalent to the 
original medium, according to a determined type of physical transfer function (known for 
the transposed medium). 

Please replace paragraph [0010] as follows: 

[001 0] For any elementary volume of the reservoir, the pressure of the oil contained 
in this volume is governed by the following equation, if gravity is not 
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taken into account : 
where: 
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§ : represents the pore volume, 

C T , the total compressibility (fluid + rock), 

K, the permeability of the rock, 

|j, the viscosity of the fluid, 

Q, the incoming flow, which is zero everywhere except in places where the well 
communicates with the reservoir, and 
P, the unknown pressure. 

Please replace paragraph [0020] as follows: 
DETAILED DESCRIPTION OF THE INVENTION 

1) Fracture network discretization 

[0020] The method comprises meshing the multilayer fractured reservoir modelled for 
example by means of the fractured reservoir modelling technique described notably in 
the aforementioned French Patent 2,575,9 4 7 2,757,947 assuming that the fractures are 
substantially perpendicular to the layer boundaries. 

Please replace paragraph [0026] as follows: 

a) Block geometry 

[0026] In order to determine the geometry of the blocks in a given layer, a two- 
dimensional problem has to solved which finds, for each fracture mesh, the points that 
are closer to this mesh than to another one . This problem is solved by using for 
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example, but preferably, the geometric method described in the other aforementioned 
French Patent 2,757,9 4 7 2,757,945 , which allows, by discretizing the fractured layer into 
a set of pixels and by applying an image processing algorithm, to determine the 
distance from each pixel to the closest fracture. During the initialization phase of this 
algorithm, value 0 (zero distance) is assigned to the pixels belonging to a fracture and a 
high value is assigned to the others. Furthermore, if the number of the fracture mesh to 
which each fracture pixel belongs is given in this phase, the same algorithm finally 
allows to determine, for each pixel : 

- the distance from this pixel to the closest fracture mesh, 
the number of the closest fracture mesh. 

Please replace paragraph [0030] as follows: 
b) Matrix-fracture transmissivity 

[0030] On the assumption of a pseudosteady state where the exchange flow is 
considered to be proportional to the difference in the pressures of each one of them, 
the following relation exists: 




where : 



F mf is the matrix-fracture flow, 



T mfl the matrix-fracture transmissivity, 



(j, the viscosity, 
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P m , the pressure of the matrix block, and 

Pf, the pressure of the associated fracture mesh. 

Please replace paragraph [0032] thru [0034 as follows: 

[0032] For this calculation, it is assumed that, in the matrix block, the pressure varies 
linearly as a function of the distance from the point considered to the fracture mesh 
associated with the block. The transmissivity is defined as follows : 



where : 

I is the length of the fractures in the fracture mesh, 
H, the height of the layer (and of the matrix block), 
K, the permeability of the matrix, and 

D, the distance from the fractures for which the pressure is the mean pressure of the 
block. 

[0033] I, H and K are known (factor 2 comes from the fact that the fracture has two 
faces). Calculation of D is made from the information provided by the image processing 
algorithm concerning the distances from the pixels to the closest fractures. In fact, 
according to the hypothesis of linear variation of the pressure as a function of the 
distance to the fracture, the relationship is : 
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where S is, in 2D, the surface area of the matrix block and d(s) is the function 
giving the distances between the points of the matrix block and the associated 
fracture mesh. 

[0034] In terms of pixels, the relationship is : 

where N is the number of pixels of the matrix block and d n the distance from pixel n to 
the fracture mesh. 

Please replace paragraph [0040] as follows: 

[0040] Passage from the time n (where all the pressure values are known) to the 
time n+1 is achieved by solving, in all the meshes and blocks of the reservoir, the 
following equation which corresponds to relation 1 once discretized : 

where : 

i, the number of the mesh or of the block, 

j, the number of the mesh or of the block next to i, 

P n , the pressure at the time n, 

t n , the time at the time n, 

Qj, the flow entering i, and 

Tjj, the connection transmissivity between i and j. 
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Please replace paragraph [0044] as follows: 
[0044] In fact, for a matrix block, equation (7) is expressed as follows : 

</> m c Tm + ^ </r 1 - p; +l ) - o (k 

subscripts m and f denoting the matrix and the fracture respectively, and C Tm being the 
total compressibility of the matrix. 
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